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Abstract. Based on a general unitary operation of two-qubit system, we investigate 
the dynamics of entanglement, entropy and purity. Distinctive features in the different 
initial states settings of an interacting bipartite system are explored. The striking 
effects of initial product states on entanglement shows that zeros entanglement can be 
obtained during the development of the interaction. 
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1. Introduction 

To test fundamental quantum concepts and implement various potential applications, 
including sensitive detection and quantum information processing, the engineering of 
quantum states has attracted considerable attention in recent years [lj. On the other 
hand, entanglement takes an important part in manipulating with information, storing, 
as in quantum coding [21 [20], sending, as in teleportation [3]. So, generating and 
quantifying the entanglement is one of the most important tasks in quantum information 
context. In quantum computer, qubits are the arithmetical units [I], which makes a new 
type of dynamics able to solve different problems faster than classical computer. It is 
desirable then to discuss the behavior of these qubits when a switch on/of the computer 
is taken place. 

Investigating the dynamics of the purity of quantum system is one of the most 
important topics in the context of quantum information theory [51 16] . It can be used as 
a measure of the degree of entanglement for mixed state (7[ E]. Bagan et. al. have been 
estimated the purity of a large number N of copies of a qubit state |9J. Our aim in the 
present article is twofold: The first is to discuss the dynamics of two separable qubits 
under the effect of the unitary operation. The second is to treat the phenomenon of 
swapping information from a sender (Alice) to a receiver (Bob). We suggest that Alice 
has a completely polarized qubit in x-direction, while Bob's qubit takes any form. We 
show that the purity of one qubit can purified on the expanse of the other one. 

The rest of the paper is arranged as follows: Section 2, is devoted to describe 
the general two-qubit pairs by means of its Bloch vectors and a cross dyadic, which is 
defined by 3 (X 1 3 matrix. In Section 3, we obtain the density operator for the individual 
subsystems, which can be used to read out how much information lost or gained. Several 
examples are given in Section 4, in which the initial states are assumed to be product 
states. In section 5, we consider a class for entangled states. The paper ends with a 
conclusion in Sec. 6. 

2. Two qubit pairs 

Assume that the user Alice and Bob are given a state of a two qubit pairs to perform 
any quantum information task. The general form of the 2-qubit state is 

p ab = \{l + 8 ■ a 1 + t ■ T l + a ■ i ~C ■ t 1 ), (1) 

where a and r are the Pauli's spin vector of the first and the second qubits respectively 
[TO] . The statistical operator for the individual quibits are specified by their Bloch 
vectors, s = (?) and t = (r). The cross dyadic C is represented by a 3 x 3 matrix, it 
describes the correlation between the first qubit, p a = tr b {p ab } = |(1 + s ■ a^) and the 
second qubit p b = tr a {p ab } = |(1 + t ■ r ). The Bloch vectors and the cross dyadic are 
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Now, we suppose that the interacting Hamiltonian of the two qubits is given by 

3 3 3 

8=1 j"=l i,j=l 

The first two terms represent the free evolution of the individual two qubits, while the 
third term represents the evolution of the correlation part. It is described by a dyadic 
consists of 9 elements. By using a suitable local unitary transformation, we can always 
choose a reference system in which this dyadic can be reduced to a diagonal one. In 
this case we have only 3 elements. In our treatment we shall consider the effective 
Hamiltonian only. In the computational basis the unitary operator [TT] is given by 

u = ri(|oo)(oo| + |ii)(n|) + r 2 (|oi)(oi| + |io)(io|) 
+ r 3 (|io)(oi| + |oi)(io|) + r 4 (|ii)(oo| + |oo)(ii|) (4) 



where, 
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This is canonical unitary operator [T2l [13] . Under the canonical evolution the density 
matrix of the 2-qubit state transforms to 
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where the new Bloch vectors and the dyadic are given by, 
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cos 2a3(sa; cos 2«i + c zy sin 2«i) + sin 2a 3 (i j: sin 2ai 



cos 2c*i), 



(1 — 2 cos 2 qj 3 ) (sy cos2a 2 — c«x sin2o;2) — sin2a 3 (i y sin2a 2 + c xz cos 2a 2 ), 

s z [cos 2 (ai + a 2 ) + cos 2 (ai — a 2 ) — 1] + t z [cos 2 (ai + a 2 ) — cos 2 (a\ — a 2 )} 
1 , 
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sin 2a 3 (s x sin 2^! — c zy cos 2a!) — (1 — 2 cos 2 a 3 )(t x cos2a! + C yz sin 2a±), 
— sin 2a 3 {s y sin2a 2 + c zx cos 2a 2 ) + (1 — 2 cos 2 a 3 ){t y cos 2 a 2 — c xz sin 2a 2 ), 
) — cos 2 (ai + a 2 )) + t z {cos 2 (a.i — a 2 ) + cos 2 (ai + a 2 ) — l) 
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(9) 
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Cll c xx , 

C12 = - ^(s 2 + t 2 ) [sin2(ai - a 2 ) - sin2(a 1 + a 2 )\ , 

- (c xy + c yx ) [cos 2 («i - a 2 ) + cos 2 («i + a 2 ) - l] , 

c i3 = (sj/ — Czx) sin2a 3 sm2a 2 — (1 — 2 cos 2 0:3) (t y sin2a 2 + c xz cos 2a 2 ), 

c 2 i = - ^(s 2 + t 2 ) [sin2(ai - a 2 ) - sin2(ai + a 2 )} 

- (c xy + c yx ) [cos 2 («i - a 2 ) - cos 2 («i + a 2 )] , 

C22 = Cyy 

C23 = — sin 2as(s x cos 2ai — c zy sin 2a±) — (1 — 2 cos 2 a^)(t x sin 2a>\ — c yz cos 2ai), 

C31 = — sin 2as(t y cos 2ct 2 — c xz sin 2a 2 ) — (1 — 2 cos 2 0^3) (s y sin 2a 2 + c 2X cos2a 2 ), 
C3 2 = sin 2a 3 (t E cos 2ct\ — c yz sin 2a±) — (1 — 2 cos 2 a^)^ sin2«i — c 2J/ cos2«i) 
C33 = c 2Z . (10) 

3. The dynamics of the individual qubits 

In this section, we investigate how the density matrices of the two qubits p a and pb 
evolve. Tracing out the second qubit (Bob's qubit), we obtain the first qubit (Alice's 
qubit),as 



/>« = i(i + iV), (11) 



where 



s~i = cos2tas(s x cos 2toi + c zy sin 2a±) + sin 2ta 3 (t r sin 2ta\ — c yz cos2toi), 
s 2 = cos2ta^{c zx sin2to; 2 — s y cos2ta 2 ) — sin2to 3 (t 2/ sin2a 2 + c X2 cos2ta 2 ), 

s~ 3 = ^(s 2 + £ 2 ) (cos 2i(o!i + a 2 ) + cos2t(a! - a 2 )) 

- ^{c xy + c ?/:r )(sin2t(ai + a 2 ) + sin2t(ai - a 2 ), 

(12) 

Similarly, if we trace out the first qubit we obtain the density matrix for the second 
qubit. In this context, we can write pb in the following form, 

h=\il+t -r l ), (13) 

where, 

t\ = sin 2taz(s x sin 2tct\ — c zy cos2tai) + cos 2tas(t x cos 2ta\ + c yz sin 2tai), 
t 2 = — sin 2ta 3 (s y sin 2ta 2 + c zx cos 2ta 2 ) —cos2ta 3 (t y cos2ta 2 + c xz sm2ta 2 ), 
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t 3 = -(s z + t z ) [cos 2t{oti + a 2 ) + cos 2£(a 1 — a 2 )] 

+ ^{c xy + c yx ){sm2t(ai + a 2 ) - sin2t(ai - a 2 )}, 

(14) 

Now, we have all the details to study the evolution of any physical property of the total 
system which is represented by p a b and its individual systems p a and p&. For example 
to quantify the amount of entanglement contained in the entangled states, we shall use 
a measurement introduced by K. Zyczkowski [TH [151 US]- This measure states that if 
the eigenvalues of the partial transpose are given by Xj, j = 1, 2, 3, 4, then the degree of 
entanglement, DOE is defined by 

4 

DOE = J2\ X j\~ 1 ( 15 ) 

i=j 

Also, we can study the evolution of the purity of the individual subsystems which 
is given by, 

rji = trpf, (16) 

where i = a,b (see for example [TTJ, HE]). Further, the amount of information contained 
in a density operator p is defined by the Von Neummann entropy S = — In A^, 
where Aj are the eigenvalues of the density operator [19J. If we start with an entangled 
initial state p a b, then we search for the survival amount of entanglement under the effect 
of the Hamiltonian. In what follows we shall consider two classes: The first for product 
states and the second for entangled pure state. 



4. Product Class 



The first class: In this section, we consider the evolution of the two qubits, where the 
initial two states are given by 

Pa = g(l + s x^x) and p b = -I. (17) 

Under the unitary transformation (jlj), the product state , p a ® Pb of the two qubits ffTTj) 
evolves to its final state with density operator, 

Pab = ^(1 + s x cos(2toi) cos(2to 3 )a 1 . + s x cos(2toi) sin(2ta 3 )r x ) 

— s x cos 2tot\ sin 2ta 3 a y T z — s x (l — 2 cos 2 ta 3 ) sin 2taia z r y ) (18) 

In Fig. (la), we plot the degree of entanglement as a function of the scaled time. In these 
calculations we assume that Alice and Bob perform the same local unitary operation. 
Also the unitary operator is defined where «i = a 2 = a 3 = | and s x = 1. It is clear 
that as soon as the interaction is switched on, an entangled state is generated. For any 
value of identical «j, this state swap from the entangled into separable states. But as 
«i decreases, the robust of remaining in entangled or separable state is decreased. 
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Figure 1. (a) shows the degree of entanglement, (b), shows the purity for the output 
states p a (solid- line) and p D (the dot-line), (c), shows the entropy for p a (solid- line) 
and p (the dot-line). The parameters are, ai = 5; s x — l,s v = s z = and 

t X ty t Z 0. 



Since, we are interested in the behavior of the individual qubits through the 
evolution, we calculate the density operator for each qubit. Consequently, 

1 1 
pi = -(1 + cos2taicos2ta 3 a x ), p~ 2 = -(1 + sm2tcti sm2ta 3 T x ). (19) 

For these new states of the two qubits, we investigate how the purity and entropy 
evolve. In Fig. (16), we consider the purity as a function of the scaled time, for the same 
parameters as in Fig.(l). From this figure we see how the purity is completely swapped 
from one pair to another. This phenomena is very important in quantum information 
context. It is possible to Alice to send all the information she has. We assume that Alice 
wants to distribute a quantum key to Bob for performing a quantum cryptography. She 
will code the secure key in her qubit by applying any coding protocol. Then she agrees 
with Bob on the unitary operators which they have to do. Fig.(lc), shows the behavior 
of the entropy for each qubit. It is clear that, the phenomena of swapping is satisfied 
for the entropy, where as the entropy of one qubit vanish the purity of the second qubit 
is maximum |20j . 

Second class: For this class, we assume that Alice and Bob both have the same 
qubit but on a different polarization as: 

Pa = ~(l + a x ), p b = -(l-r x ), (20) 

Namely s x = l,t x = — 1 and the other components vanish. After the evolution, the 
joined density operator is defined by the Bloch vectors 

s = (cos2t9 1 cos2t6> 3 - sin2t9 l sin2£0 3 , 0, 0), 
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t = (cos 2^i(l - 2 cos 2 te 3 ), -cos2t0 2 (l - 2cos 2 t0 3 ), 0), (21) 
in addition to the dyadic^, which is defined by 

C X y C XZ 0, 

c yx = c yy = 0, c yz = sin 2*0^1 — 2cos 2 2t9 3 ) — cos2t9 1 sin2£0 3 , 

c zx = 0, c zy = - sin 2t0i (1 - 2cos 2 2t9 3 ) - cos 2^ sin2£0 3 , c zz = 0, (22) 

Also, the individual qubits are specified by their new Bloch vectors 

s = (cos 2t0i cos 2t0 3 , 0, 0), 

t = (-cos 2t0i cos 2*0 3 + sin 2t9 1 sin 2t0 3 , , 0). (23) 

In Fig. (2), we plot the degree of entanglement against the scaled time. From this figure, 
it is clear that the generated state is entangled most of the time except at specific 
points. At these points the two qubits are completely separated. This is clear from 
Fig. (26), where we plot the purity of the two qubits. We see when the purities of 
the two qubits are completely swapped, the two qubits become separable states. This 
occurs at normalized t = 1.5,2.5,3.5 etc. On the other hand, when the purities of the 
two qubits coincide, the generated entangled state has a large degree of entanglement 
and it reaches its maximum value when the purities are maximum. This phenomena 
is seen once the interaction starts as well as at normalized t = 3n, n = 1, 2, 3, .... This 
feature is confirmed from Fig(2c) , where the entropies of the two qubits are plotted. In 
some intervals neither Alice nor Bob know any information about the other. At these 
instance all the information has been transformed from one pair to each other. Alice 
can employ this technique to send information to Bob. 
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Figure 3. The same as Fig.(l), with t x = t y = 0.5 and t z = 0. 



Third class: In this class, we assume that Alice has the same qubit as before, while 
Bob's qubit is polarized along x and y axis. In explicit form 
1 1 

Pa= ^(1 + S X (T X ), p b = -(l+t x T x +t y T y ). (24) 

After the evolution the joined density operator p ab is defined by the following Bloch 

vectors 

s = (s x cos 2t9 1 cos 2t9 3 + t x sin 2t9\ sin It sin 3 , —t y sin 2t9\ sin 3 , 0) , 
t = (—t x cos 2t9i (1 — 2 cos 2 t9 3 ) + s x cos 2t9\ sin 2t sin 9% , 

-t x cos2t9 2 (l -2cos 2 9 3 ), 0), (25) 

and the cross dyadic 1 C 

c xx = c xy = 0, c xz = -t y sin 2^ 2 (1 - 2 cos 2 £0 3 ), 

= c ra = 0, c yz = — t x sin 2i0i(l — 2 cos 2 10 3 ) — s x cos 2£0i sin 2t0 3 , 
c zx = — t y cos 2t92sin2t9 3 , 

c zy = - s x sin 2*0i (1 - 2 cos 2 t9 3 ) + * x sin 2t0 3 cos 2t9 u c zz = 0. (26) 
On the other hand after the evolution, the two qubits evolve as 

p a = -(1 + s x cos 2t9 x cos 2t9 3 a x — t y sin2t9 3 a y ), 

p b = i(l + (t x cos2t9 1 cos 2t0 3 + s x sin 2t9\ •sm.2t9 3 )T x 

- t y cos 2£0 2 cos 2t9 3 r y ), (27) 

In Figs. 3a, we have plotted the amount of entanglement contained in the generated 
state p a b- We consider s x = 1 while t x = t y = 0.5, and assume that Alice and Bob 
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both perform the same unitary operator with |. We can show that the amount of 
entanglement is larger than the entanglement which is depicted in the first class, but 
less than the second class. Also the purity of the initial state of Bob is larger than the 
first class and smaller than the second one. In Fig. 3c, the behavior of the entropies of the 
individual qubit is shown. This figure depicts a partially swapping of the purity. This 
is clear where the purity of Alice's qubit never reaches zero. It decreases with amount 
enough for Bob's qubit to reach the maximum purity. The amount of information which 
is known by Alice and Bob about the qubit of each other is plotted in Fig. 3c. It is clear 
that in this class, Bob is not able to know all the information contained in Alice's qubit. 

5. Entangled Class 

We consider a class of pure entangled states defined by its joint density operator 

Pab = +po- x - pr x - c xx - q(c yy + c zz )), (28) 

where q = a/ 1 — p 2 . This class of entangled state is very important, where for p — 0, one 
get the singlet state and separable for p = 1. This class has been studied in [2T] . where 
we quantify the amount of exchange information between this state and environment. 
Under the effect of the unitary operator, this state is transformed such that, 

s = (pcos2t(0 1 + 3 ),O,O), 

t = (p cos 20! ( 1 - 2 cos 2 2t9 3 ) + p cos 2t9 1 sin 2t9 3 , -p cos 2t9 2 (1-2 cos 2 t9 3 ) 

and 

Cxx 1) C-xy Cxz 0) 

Cyx — 0, Cyy — — 1, c^ = psin2t6> 1 (l — 2 cos 2 £0 3 ) — pcos2t0 1 sin2t6> 3 , 

c zx = 0, c zy = -pcos2t(6> 1 + 9 3 ), c zz = -q- (30) 

On the other hand, if we trace out one qubit, we can obtain the other, so one can obtain 
the density operators for Alice and Bob qubis as 

Pa = t;(1 +pcos2t9 1 cos2t9 3 a x ), p b = ^(1 - p cos 2t(9 x - 9 3 )t x ). (31) 

The behavior of the amount of entanglement contained in the entangled state, which 
is defined by f T2"9"|) and fl3"U|) . is shown in Fig(4a). We consider that the initial state has 
p = 0.7 as an initial state, where this parameter completely defines the state. First of all 
the state is always entangled and the degree of entanglement increases and decrease but 
never reach zero. In Fig.46 and 4c, we plot the behavior of the purities and entropies 
of the individual systems, p a and pb, respectively. From these two figures, we can see 
that the purity of one qubit may be increased on the expense of the other qubit. On 
the other hand the entropy of each of them can not reach zero. This means that, the 
two qubits still always have some information about each other. 
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Figure 4. The figure is evaluated for p = 0.7 and q = y/l — p 2 . (a) the degree of 
entanglement for p a b, after the evolution (b) The purities for the Alice 's state p a { the 
solid line) and Bob'state p& (dot -line) . (c) The Entropy for p a (solid line) and pb (dot 
line) 



6. Conclusion 

In this contribution, we study the interaction of a two separable qubits by means of 
the canonical unitary operator. An analytical expression of the joint density operator 
is obtained. We calculate explicitly the reduced density operator for each qubit. 
Different classes of initial density operators are considered. The amount of entanglement 
contained in the joint entangled state is quantified by using the negativity measurement. 
The swapping phenomena of purity and entropy is investigated for the reduced density 
operator. It is clear that the unitary operator parameters ccj, i = 1,2,3 play the central 
role in the interaction process. If the state of one qubit is different from the other qubit 
or polarized in an opposite direction, then the information is completely transformed 
from one qubit to the other. 
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